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Abstract 

We consider the neutralino mass matrix in a supersymmetric theory based 
on SU(2) X U(l) gauge group with arbitrary number of singlet, doublet and 
triplet Higgs superhelds. We derive an upper bound on the mass of the 
lightest neutralino, and a lower bound on the mass of the heaviest neutralino, 
in such a theory. Assuming grand unihcation of the gauge couplings, the 
upper bound on the mass of the lightest neutralino can be expressed in terms 
of the gluino mass. For a gluino mass of 200 GeV, the tree level upper bound 
on the mass of the lightest neutralino is 92 GeV, which increases to 166.5 
GeV for a 1 TeV gluino. We also discuss the effect of dominant one-loop 
radiative corrections on these bounds. 
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In supersymmetric gauge theories, each fermion and boson of the Stan¬ 
dard Model is accompanied by its supersymmetric partner, transforming in 
an identical manner under the gauge group []^. In supersymmetric theories 
with R-parity conservation [Q], the lightest supersymmetric particle (LSP) is 
expected to be the lightest neutralino, which is the lightest mixture of the 
fermionic partners of the neutral Higgs and neutral electroweak gauge bosons. 
In order to give masses to quarks and leptons, and to cancel triangle gauge 
anomalies, at least two Higgs doublets Hi = and H 2 = {H^jH^), 

with opposite hypercharge iY{Hi) = — 1 ,Y{H 2 ) = +1), are required in the 
minimal version of the Supersymmetric Standard Model (MSSM) [|^. The 
fermionic partners of these Higgs bosons mix with the fermionic partners of 
the gauge bosons to produce four neutralino states i = 1,2, 3,4, and two 
chargino states yf, i = 1,2, in the MSSM. In the nonminimal supersymmet¬ 
ric standard model containing a Higgs singlet, besides the two Higgs doublets 
of the minimal model, the mixing of fermionic partners of neutral Higgs and 
gauge bosons produces hve neutralino states. The neutralino states of the 
minimal model [^, ^ and those of the nonminimal model isa have 

been studied in great detail, because the lightest neutralino, being the LSP, 
is the end product of any process involving supersymmetric particles in the 
hnal state. 

Recently we considered the lightest neutralino state in a general super- 
symmetric theory containing an arbitrary number of singlet, doublet and 
triplet Higgs superhelds under the Standard Model gauge group, and ob¬ 
tained an upper bound on its mass p[|. If we assume the simplest form of 
grand unihcation of gauge couplings in such a theory, whereby the triplet 
and extra doublet Higgs helds are eliminated, then this upper bound is con¬ 
trolled by a soft supersymmetry gaugino mass (which we can take to be the 
gluino mass), and the vacuum expectation value of the doublet, but not the 
singlet, Higgs helds. Since the latter are known, the bounds are controlled 
by the gluino mass alone. For a gluino mass of 200 GeV, the upper bound 
on the lightest neutralino mass (including the dominant one-loop radiative 
corrections) was shown to be 62 GeV, which increases to 178 GeV for a 1 
TeV gluino. 

In this paper we generalize the above result. We consider the neutralino 
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mass matrix in a general supersymmetric theory with an arbitrary Higgs sec¬ 
tor, which includes singlets, doublets and triplets under the Standard Model 
gauge group. We obtain an upper bound on the mass of the lightest neu- 
tralino and a lower bound on the mass of the heaviest neutralino state in 
such a general supersymmetric theory. These bounds depend on the soft 
supersymmetry breaking gaugino masses and the vacuum expectation val¬ 
ues of the doublet and triplet, but not the singlet, Higgs helds. This is in 
contrast to the situation that obtains in the Higgs sector of such a general 
supersymmetric theory |P], where the (tree level) upper bound on the light¬ 
est Higgs boson mass is controlled by the vacuum expectation value of the 
doublet Higgs helds and dimensionless parameters only. Nevertheless, if we 
assume grand unihcation of such a general supersymmetric theory, then we 
show that the upper bound is controlled by a soft supersymmetry breaking 
gaugino mass parameter (which can be taken to be the gluino mass), and 
Mz and 9w Since the latter are known, the former entirely determines the 
upper bound. 


We start by recalling the neutralino mass matrix in the Minimal Super- 
symmetric Standard Model |]l|. In the basis 

V’J = {-ix', j = 1, 2, 3, 4, (1) 

where X! and are the two-component gaugino states corresponding to the 
U(l)y and the third component of SU(2)i gauge groups, respectively, and 
'ipfjl are the two-component Higgsino states with hypercharge -1 and 

-|-1, respectively, the neutralino mass matrix can be written as 


M = 


Ml 0 

0 M 2 

—Mz sin^vi^ cos/3 Mz cos9w cos(3 
Mz sin sin /3 — Mz cos 9w sin P 


—Mz sin^vr cos/3 Mz sinsin/3 
Mz cos 0vrcos/3 — Mz cos 9 w sin P 

0 —fj, 

—g 0 


( 2 ) 


where Mi and M 2 are supersymmetry breaking gaugino masses associated 
with the U(l)y and SU(2 )l subgroups of the standard model, respectively, 
and fi is the Higgs(ino) mixing parameter in the superpotential. Here, 
tan/3 = V 2 /V 1 , where vi = {Hf) and V 2 = ^^e vacuum expec¬ 

tation values of the neutral components of the two Higgs doublets, and 


3 




— ( 5 '^+ 5 '^^)('^i+'^ 2 )/ 2 , with g and g' being the gange conplings associated 
with the SU(2) and U(l) factors of the standard model gange gronp, respec¬ 
tively. Diagonalization of the mass matrix (|^) gives the physical nentralino 
masses and eigenstates [|T^]. We shall denote these nentralino eigenstates by 
Xi, X 2 ) Xs) cind X 4 , labelled in order of increasing mass. Since some of the 
nentralino masses resnlting from diagonalization of the mass matrix can 
be negative, it is convenient to consider the sqnared mass matrix MlM and 
the corresponding sqnared masses of the nentralinos. An npper bonnd on the 
sqnared mass of the lightest nentralino Xi can be obtained by nsing the fact 
that the smallest eigenvalue of is smaller than the smallest eigenvalne 

of its npper left 2 x 2 snbmatrix 


Ml + Ml; sin^ 9w — Af| sin 9w cos 9] 
—M| singly cos 6 'vi/ M|-|-M| cos^ 


'w 


(3) 


thereby resnlting in the npper bonnd |Tl 


< min(M 4 -|- M| sin^ 9w, M^ -|- M| cos^ 9w) 


(4) 


On the other hand, the larger eigenvalne of the npper left 2x2 snbmatrix (|^) 
of iVflM gives a lower bonnd on the sqnared mass of the heaviest nentralino: 

> max(M 4 -|- M| sin^ 9w, Af| -|- M| cos^ 9iy). (5) 


We note that the bonnds (^) and (j^) are controlled by, in addition to Mz 
and 9w, the soft SUSY breaking gangino masses. Mi and M 2 . This is in 
contrast to the Higgs sector of MSSM, where the corresponding bonnds on 
the (tree level) masses of the lightest and the heaviest scalar Higgs bosons 
are controlled by Mz, and not by snpersymmetry breaking masses H- 


We now consider a general class of supersymmetric standard models, 
namely, supersymmetric models based on Standard Model gauge group with 
an arbitrary Higgs sector. We shall assume 
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(i) Two Higgs doublets with hypercharge Y = =pl, which are 

coupled to the quarks and leptons in the superpotential 

1T° = huQLUlHi^^ + hnQLDlHi^^ + hElLElH^^^. (6) 


In addition, we assume an arbitrary number of extra pairs of Higgs 
doublets Hi\H 2 \j = 2, 1, which are decoupled from quarks 

and leptons so that there are no dangerous flavor changing neutral 
currents 0. 


(ii) Gauge singlets = 1, ...,ns. 

(hi) SU(2) triplets with Y = 0. 

(iv) SU(2) triplets with Y = ±2. 


Since the two Higgs doublets and are the minimum number of dou¬ 
blets which are required in MSSM to give masses to all fermions and to cancel 
triangle gauge anomalies, the above extra Higgs multiplets are the only ones 
that can provide renormalizable couplings to all possible combinations of 
and H 2 ^\ viz. H 2 ^\ h[^\ and H 2 ^ H 2 \ in the superpotential. 

More exotic Higgs representations contribute to the /d-functions of the gauge 
couplings and give lower values of the upper bound on the lightest Higgs 
mass in general supersymmetric theories. Since we are interested in absolute 
upper bounds on the particle masses in these theories, we do not consider 
such representations. The most general renormalizable superpotential for the 
above Higgs supermultiplets can be written as 


W' = Wi + W 2 , 


(7a) 


where (repeated indices summed) 
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halm ^ 

6 

, prill) pT{u) pT(a) 

6 


(7c) 


and where we have represented the triplets by 2 x 2 traceless complex ma¬ 
trices, and the trace in is over the matrix indices. More explicitly, the 
triplets are represented as (with mnltiplicity indices suppressed): 


e/V2 

r -e/V2 


(8a) 


■ ++/ v ^ 


, T2 = 

■ ++/+2 



-++/+2 _ 

^2 

-++/ V 2 _ 


(8b) 


Without loss of generality, we can make a unitary transformation in the 
space of Higgs doublets and such that only the Higgs doublets 


and aquire |T^ a non-zero vacuum expectation value ((if)^'*'^) = Vi, 
(772^^°) = ^ 2 )- This requires that some of the Yukawa couplings in (|7b|) 
vanish: 


( 1 ) 0 \ _ 


=/i’“=st”=sj'”= 0 u^i). (9) 


We will assume this condition in what follows. With this condition, the part 
of the superpotential (|7aD involving only the neutral helds can be written as 


W' = 


-E 

(7 

+ E 

a 

-E 


d+l 

pier pjii)0 j^{i)0 

i=l i^i,i¥=j 

d+l 

^/+77{*^°77f^ 

i=i 

d+l 

^ (7+77P°77®° + ^ +™77f°77^^')° 

i = l ir/rjr/rl 




^(a)0 

Tf 




(m)0 


6 

























d +1 


2 = 1 


-E 

m 

+E 

a 

+E 

{il'cr 

Choosing the neutral fields 


^(a )0 




(m)0 


. Im 

Ehfl _/Y (C _/Y (C _/Y O 
6 




( 10 ) 


'4^1 ( 1 JjU)0 ^ JJ^)0 , Ip jj(j)0 , ijj ljj^(a)0 , 'ljj^(l)0 , ijj^Wo'j , 

j = 2, ....,d+1, a = l,....,ns, a = l,....,to, 1 = 1, ...,ti, ( 11 ) 

/ = 1,., (4 + 2d + + to 4“ 2ti), 

as the basis, we can write the neutralino mass matrix for the general class 
of supersymmetric models based on standard model gauge group with an 
arbitrary Higgs sector as 


M 


Mi 3 

0 

Mi 4 

0 

0 

0 

2V2g'y[ 

1 

to 

_1 

M23 

0 

M24 

0 

0 

gV2u‘^ 

-gV2y\ 

gV^yi 

2g[y[ 

0 

f 

0 

/r«2 

f^V2 

29(^1 

0 


1 

to 

P 

fij 

rpv2 

fpV2 

0 

0 



2922/2 

0 

/r«i 

f^vi 

0 

292^/2 

SYMM 



-292^*92 

0 

0 

0 

0 






0 

0 

0 






0 

Ki^vT 

h' ,yT 

aml^ 1 







0 

h' ,n“ 

ami 


0 


( 12 ) 

where 

Mi 3 = — M^sin 6 *vi/cos/5, M 14 = sin sin/?, 

M 23 = Mz cos 9w cos P, M 24 = — Mz cos 6w sinP, 
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are the matrix elements which enter the nentralino mass matrix (§) of the 
MSSM, and we have defined the various vacuum expectation values and the 
couplings as follows: 

«“ = ({<“>"), x'’={/V<''>>, = (13a) 

/r = / 2 “ = /“ 7 A /r =/f = /2’7A (ub) 


/ = ifix” + />“). /■ = + /“«“). /« = + /«“ u 772 ), 

(13c) 

a[ = 92 = - 9 “'. Kim = (13d) 


In (|13aj) we have represented the Higgs components of the superhelds by the 
same symbol as the superhelds themselves. The upper bound on the squared 
mass of the lightest nentralino can be obtained by examining the upper left 
2x2 submatrix of M corresponding to M of This submatrix can 

be written as 


■ + M| sin^ 0VC + 6c/'V 

—M| sin 6w cos 6w — “^gg'u^ 


—M| sin 6w cos 6w — ‘^gg'y^ 
M| + M| cos^ 6w + ‘2g‘^u‘^ 


(14a) 


where 


V- = EKsi)" + (A)% K = (14b) 

l a 

and where we have used the expressions for W and Z masses appropriate for 
the general theory that we are considering: 

= 2 {g‘^ ^g'‘^)[^i'^ 2 (iSa) 
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Mw = + V]- 


(15b) 


The vaccum expectation values that enter into the W and Z masses 
and (|15b|) are experimentally constrained by the p parameter. From a recent 
global fit, which includes the CDF data, we have m 


p= 


Ml 


Ml COS^ 9]y 


= 1.0002 ±0.0013 ±0.0018, 


(16) 


where the second error is due to the Higgs mass. This result is remarkably 
close to the expected Standard Model value of p = 1. Taking a value of p = 1 
implies, through (|15a|) and (|15b|) , the following relation between the triplet 
vacuum expectation values: 

= 2y^. (17) 

We shall henceforth assume (0 to be true. With constraint (O). the W 
and Z masses can be written as 0 

A/J, = Af| cos“ + I'l + 

Thus, the combination (vf -hv^ ±4p^) is determined to be ~ (174GeH)^, but 
the ratio of the triplet to the doublet vacuum expectation values y/ (n^±n|)^/^ 
is not fixed. Using (0), the 2x2 submatrix ([14a|) can be written as 


Ml ± M| sin^ 6w ± ~Mz cos 9w sin 9w — ‘^gg'y^ 

-M| cos 9w sin 9w — ‘^gg'y^ M^ ± M| cos^ 9w ± g^y"^ 


(18) 


The smallest eigenvalue of the 2x2 matrix ([Tq) serves as the upper bound 
on the squared mass of the lightest neutralino in the general supersymmetric 
Standard Model: 


M^q^ < min(M^ ± M| sin^ 9w ± M| ± M| cos^ 9w ± g‘^y‘^)- (19) 

On the other hand, the heaviest neutralino (Xn) mass is bounded from below 
by 

M^o > max(M^^ ± M| sin^ 9^ + 6g'‘^y‘^, M| ± M| cos^ 9w ± g^yl- (20) 
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These bounds depend on a priori unknown vacuum expectation values of 
the triplet Higgs fields i '^ 2 '^) > ^ind the supersymmetry breaking gaugino 
mass parameters. Nevertheless, as we shall see these bounds can become 
meaningful in theories with gauge coupling unification [ITl]. We note that 


the singlet vacuum expectation values decouple from these bounds. 


In the general supersymmetric theory that we are considering, the renor¬ 
malization group equations (RGBs) [|^ for the standard SU(3) x SU(2) x 
U(l) gauge couplings can be written as {gf = g^ = g^, tanOw = 

9 '19) 93 is the SU(3) gauge coupling constant) 


Ibvr 

Ibvr 


2 ^9i 
dt 

2^92 


33 

5 


3 

5 ' 


— + 7 ( 6^1 + d) gf, 


dt 


— [1 + 2 to + -|- d)] g2, 


2^93 


Ibvr —— = 
dt 


M- 


( 21 ) 


Obviously, these RGBs depend on the number and the type of the Higgs 
representations (d, toGi)- We note that the additional doublets increase the 
jd functions of the gauge couplings, even though the VBVs of the additional 
doublets have been rotated away. If we assume that the gauge couplings 
91,92, and gs unify at some grand unification scale, i.e. gi{Mu) = g 2 {Mu) = 
f 2 ,{Mjj), with Mjj the unification scale, then the simplest choice is ^ 


d = to = H = 0, 


( 22 ) 


with rig arbitrary, i.e. MSSM with an arbitrary number of singlet superhelds. 
In other words, only Higgs singlets, besides the two Higgs doublets of MSSM, 
are consistent with the simplest form of unihcation JI^.The bounds (|I^) and 
( pOD now reduce to 


^^0 < min(M^^ -|- sin^ 6w, -|- M| cos^ 6\y) 


(23) 


M^o > max(M^^ -|- sin^ 9w, -|- M| cos^ 9w), (24) 

which are identical to the corresponding bounds (^ and (^ in the minimal 
supersymmetric standard model [^. 
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However, gauge coupling unification can be achieved even when d, to, ti > 
0 by adding colored multiplets or introducing an intermediate scale in the 
theory |^. Therefore, it is important to ask what are the bounds on the 
mass of the lightest neutralino in the most general case without assuming 
any particular form of grand unihcation with a specihc particle content. To 
evaluate the upper bound (|T^) in the most general case, we recall that the 
gaugino mass parameters satisfy the renormalization group equations (RGBs) 


21 


167r" 


,dMi 

dt 


= hMisi, 


(25) 


where the coefficients bi are the same that occur in the evolution (^l|) of 
the corresponding gauge couplings. Equations (|2l|) and p5| ) imply (ctj = 
g"flA'K,au = 


Mi{Mz) / oii{Mz) — M2{Mz) / ci2{Mz) — M^{Mz) / ot'j,{Mz) — vfi\i2loi\j, 

(26) 

where M 1/2 is the common gaugino mass at the grand unihcation scale, and 
au is the unihed coupling constant. It is important to note that (p 6 D, which 
is a consequence of one-loop renormalization group equations, is valid in any 
grand unihed theory irrespective of the particle content. Two-loop ehects to 

0 


are expected to be numerically small 


relation can be violated only at order au/n [^. The relation 


We also note that the GUT 
]) is the 


same which occurs in the MSSM with grand unihcation. It reduces the three 
gaugino mass parameters to one, and since the the gluino mass rrig is equal to 
I M 3 1, it is, therefore, convenient to let the gluino mass be the free parameter. 
The other two gaugino mass parameters are then determined through 


M 2 {Mz) = {a/9w)mg — 0.27mg, 

Mi{Mz) = {5a/3aocos^ 6w)mg — 0.14mg, 


(27) 


where a is the hne structure constant, and where we have used pj, ^ 
the values of couplings at the energy as 


a 


. -1 


{Mz) = 127.9, sin^^vv = 0.2315, 03 = 0.123. 


(28) 


Using (pTl) and the most conservative estimate of the upper bound on the 
triplet vacuum expectation value, y‘^ < (87GeU)^, following from {v\ + V 2 + 
Ay'^) ~ (174GeU)^, the upper bound ( 0 ) can be written as 
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Mj; < M? + M| siii^ + 6gV ~ (0.02m| + 7733.5)GeV^ 


(29) 


Note that the lower bound (|20|) on the mass of the heaviest neutralino cannot 
be evaluated in the most general case that we are considering. For a gluino 
mass of 200 GeV, the upper bound (^) on the mass of the lightest neutralino 
is 92 GeV. Similarly, for a gluino mass of 1 TeV, the upper bound becomes 
166.5 GeV. This should be compared with the upper bound of 52 GeV and 
148 GeV, respectively, for gluino mass of 200 GeV and 1 TeV, which one 
obtains when one assumes H a simple form of gauge coupling unification 


O- 


The bounds (P!P|), (2D) and (|2D|) are our main results. While the former 
are valid in a supersymmetric theory based on a standard model gauge group 
and an arbitrary Higgs sector, the latter is valid only when we impose the 
constraint of grand unification on gauge couplings in such a supersymmetric 
theory. Thus, although a bound on the lightest neutralino mass exists in an 
arbitrary supersymmetric theory, it is calculable in terms of gluino mass only 
under the assumption of gauge coupling unification. We further note that 
since the soft supersymmetry breaking gaugino masses satisfy the GUT re¬ 
lation (p6D in any grand unified theory based on a simple group independent 
of the breaking pattern to the Standard Model gauge group iQ, our result 
is valid in any such supersymmetric grand unihed theory. 


Finally, we discuss the effect of radiative corrections on the tree level 
upper bound on the lightest neutralino mass derived above. The dominant 
radiative corrections arise from top-stop loops. If the neutralino is predom¬ 
inantly Higgsino-like, then a simple estimate of the radiative corrections to 
its mass arising from a top-stop loop is given by 


AM^o 

M^o 


~ 3 



'fnf 1 

~ 5%, 


(30) 


where the factor 3 comes from color, comes from the top Yukawa couplings 
at the two vertices in the loop (ht = see (||)), and l/167r^ is the usual 
loop factor. In making the above estimate, we have taken the squark mass to 
be equal to 1 TeV. A similar estimate holds for a gaugino-like neutralino, and 
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also yields an estimate for the correction of order 5%. These estimates are 
very close to the generic resnlts which emerge from detailed calcnlations [29 


carried ont in the context of the minimal snpersymmetric standard model, 
althongh, for some extreme valnes of parameters, the radiative corrections 
to the lightest nentralino mass can be as large as 20%. Taking these resnlts 
as indicative of the radiative corrections in the general model that we are 
considering, we estimate a conservative radiatively corrected npper bonnd 
on the mass of the lightest nentralino to be abont 110 GeV for a 200 GeV 


glnino, which increases to 200 GeV for a glnino mass of 1 TeV ||3b 
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